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FINAL TECHNICAL REPORT 


A METHOD FOR EXPERIMENTAL MODAL SEPARATION* 

I. Introduction 

A systen of cornHiter software has been developed for the numerical simula- 
tion of multiple-shaker modal survey testing. A fairly simple analytical method 
using simulated experimental data is employed to optimize the shaker force- 
anqilitude distribution for the purpose of isolating individual modes of vibra- 
tion. This method is referred to here as the RFAD method, for relative force- 
amplitude distribution. 

The system has been employed in a thorough simulation of modal testing on 
a small, fairly simple, analytical structural model. The detailed description 
of this simulation is included in a Master's degree thesis by Stafford (1976). 
The software system is currently being used to test the RFAD method on various 
other structural models, both discrete and continuous. Each model was designed 
to have at least two closely spaced modes in order that we might simulate the 
separation of close modes by multiple-shaker excitation. Inertia, damping, 
stiffness, and modal data for these models are stored on magnetic disks, 
available by direct access to the interactive Fortran programs which perform 
all computations required by the RFAD method. 

The support for this research by NASA has been renewed for an additional 
year under the same number. Grant No. NS6 1276. Hence, this report is not 
the final presentation of theory, procedures and results to be submitted on 
our research into the RFAD method. 


* The NASA Technical Officer for this grant is Mr. Robert Miserentino, 
NASA Langley Research Center. 
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II. Short Description of the RFAD Method 

The* RrAD method is based on the following theoretical result for an n- 
degree-of- freedom damped structure. In order that the structure be forced 
to vibrate in its rth undamped normal mode, the condition 

[C(o)^) ]F = 0 

must be satisfied (Asher, 1958; Bishop and Gladwell, 1963; Craig and Su, 1974). 

F is the n x 1 column vector of shaker force-amplitudes, and [C{d)^)] is the 

* 

n X n coincident-response matrix evaluated at the rth resonance frequency, 

In practice, it is generally not possible to provide forcing excitation at 
all the degrees of freedom of a structure, so the problem which the RFAD method 
considers is 

[C*]F* = 0 

where F* is a p x 1 matrix of force-amplitudes corresponding to the p available 
shaker locations, and [C*] is the incomplete p x p coincident-response matrix 
for those same locations. This is an eigenvalue problem for which, ideally, the 
eigenvalues (those frequencies at which Det [C*] = 0) are the undamped natural 

frequencies, and the eigenvectors are the force-amplitude distributions which 
will tune reasonably pure modes. 

III. Analytical Models with Closely Spaced Modes 

Perhaps the most difficult task in modal survey testing is the separation of 


* The coincident-response matrix is the real or in-phase portion of the 
complete frequency-response transfer function matrix. 
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closely spaced nodes. Hence, one of our primary objectives In developing 
analytical models to test with the RFAD method has been to simulate closely 
spaced modes. Although such modes appear frequently In vehicle modal testing 
programs, they are actually quite difficult to produce In relatively simple 
structural models of the type appropriate for this research. Hence, we have 
developed an analytical procedure for ''optimizing" a structural model to have 
at least one pair of closely spaced modes. 

The structural model we have considered most extensively is the five-degree- 
of-freedom cantilevered plane grid shown in Fig. 1. The beams are aluminum 
rods having 0.254 m (10.0 Inch) length* and 0.00635 m (0.25 inch) diameter; 



they are considered massless. All mass is considered to be concent’-ated at the 
beam intersections. Our procedure has been to adjust the concentrated masses 
in order to force together selected adjacent modes. 


* All calculations have been made In the primary units of pounds, inches 
and seconds. 
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Stafford (1976) used an early form of our analytical mode-merging method 
to produce closely spaced third and fourth modes of this ir»del having fre- 
quencies of 30.54 rad/sec and 31.90 rad/sec, respectively. However, this 
situation resulted from a set of masses weighing from 11.9 N (2.68 lb.) to 
114 N (25.7 lb.). In view of the small size of the model, the latter value 
is unrealistically large. But we try to retain a reasonable degree of 
physical realism in our analytical models, since we may elect later to build 
laboratory specimens based on them. Hence, we have refined our mode-merging 
procedure so as to hold the masses to physically realistic values. The theory 
and procedure are presented in Appendix A. With this procedure we have pro- 
duced an eigensolution for the five-degree-of-freedom model having third and 
fourth resonances of 53.50 rad/sec and 55.41 rad/sec, respectively. The masses 
in this model weigh from 3.38 N (0.760 lb.) to 20.6 N (4.62 lb.). 

IV. Simulation of Arbitrary Hysteretic Damping 

We have studied primarily n-degree-of- freedom damped models governed by 
the linear matrix equations 

[m] X + ^ [d] ^ + [k] X = f 

where [d] is the hysteretic damping matrix, and excitation and response 
are sinusoidal at frequency u, 
f = Fe'*“* 

X = Xe^“^ 

All previous studies of the RFAO method (Bishop and Gladwell, 1963; Craig 
and Su, 1974; Stafford, 1976) have considered only the case of hysteretic 
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damping which does not couple the undamped normal modes of the system, 

Stafford used a proportional damping matrix, 

[d] = g tk] 

which is uncoupled by the standard transformation to normal coordinates and 
gives structural damping constant g in all moder . 

However, there is no reason to expect in general that damping will not 
couple the modes. So we have developed the capability to simulate a structure 
with arbitrary hysteretic damping. The theory behind this simulation is 
evidently due to Bishop and Johnson (1960), and it has been discussed ex- 
tensively by Mead (1970). It involves first determination of the complex 
eigenvalues and eigenvectors of the homogeneous problem 

([k] + i [d] - X[m] ) t = 0 

and then solution for the steady-state response due to sinusoidal excitation 
by modal analysis with the complex modes. 

To provide an example of the possible effects of damping which couples the 
normal modes, we have applied the RFAO method with both proportional and non- 
proportional damping to the f1ve-degree-of -freedom model discussed in Section 
III. For the proportional case we use [d] = 0.03 [k] , and for the non- 
proportional case we use simply the diagonal portion of the proportional 
damping matrix. Since a primary objective of the RFAD method is to separate 
closely spaced modes, we seek distributions of shaker force-amplitudes which 
will separate the third and fourth modes at 53.50 and 55.41 rad/sec, respectively 
Plots of Det [C(w)] for all five masses excited and Det [C*(u))J for only 
masses 1 and 5 (see Fig. 1) excited are shown on Fig. 2. Det [C] curves 



Normalized Det [C ] Normalized Det [C] 



Frequency, (rad/sec) 

Figure 2 Coincident-Response Determinants 



for the two damping cases are barely distinguishable in the frequency range 
shown, and both have zero crossings exactly at the undamped resonances as 
is required by theory. Det [C*] curves for the two cases are similar but 
not identical anywhere, as should be expected since they are derived from 
incomplete transfer function matrices. 

Eigenvalue solutions corresponding to the zero crossings of the determinant 
plots are listed in Tables la and b. Each shaker force-amplitude eigenvector 
is normalized to its largest value. With all masses excited, there are sub- 
stantial differences between the force-amplitude eigenvectors for proportional 
and non-proportional damping. This means that the shaker force-amplitude re- 
quired to isolate a mode is strongly dependent on the specific form of damp- 
ing, as one would expect. Curiously, this dependence appears to be much 
weaker with only masses 1 and 5 excited. 

V. On the Number of Predominant Modes and the Number of Required Shakers 

One consequence of our numerical studies is considerable qualitative 
evidence of a relationship between the number of modes predominant in the 
motion at a given frequency and the number of shakers required to separate 
modes in the vicinity of that frequency. This is not a new idea, as Trail 1- 
Nash (1958), Bishop and Gladwell (1963), and Asher (1967) have previously 
discussed it. In fact, Asher proposed a quantitative technique for determin- 
ing the number of predominant modes by means of a determinant evaluation pro- 
cedure using transfer-function data. We have tested Asher's technique exten- 

* 

sively and have found it to be unreliable. 


* In contrast, Ibrahim and Mikulcik (1976) employed a similar technique with 
filtered transient response data and found it quite satisfactory. 
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a. A11 masses excited 



Proportional damping 

Non-proportional damping 

Zero crossings 
(rad/sec) 

53.50 

55.41 

53.50 

55.41 

F, 

0.324 

1.000 

0.208 

0.360 


0.070 

-0.883 

0.141 

-1.000 


1.000 

-0.896 

1.000 

-0.503 

^4 

-0.234 

-0.460 

-0.315 

-0.348 

^5 

-0.213 

0.736 

-0.060 

' 

0.116 


b. Masses 1 and 5 only excited 



Proportional damping 1 

j Non-proportional damping J 

Zero crossing 
(rad/sec) 

53.25 

55.36 

53.16 

55.31 

* 

Pi 

0.436 

0.709 

0.444 

0.719 

* 

P 5 

-1.000 

1.000 

-1.000 

1.000 


Tables la and b 


Coincident-response determinant zero-crossing frequencies and force-amplitude 
eigenvectors associated with plots in Figure 2 . 
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A concept analogous to tlie number of predominant modes is that of “a 
best approximating subspace". This has been analyzed by Cliff (1975) in 
the context of control theory. Cliff's approach appears to be applicable, 
with sonte modifications, to our problem. A detailed description of the 
theory behind this method Is presented in Appendix B. The method involves 
numerical testing of experimental or simulated experimental transfer function 
data to determine the number of predominant modes. Within the context of 
this method, there are several different possible ways to evaluate a given 
set of data. We are currently conducting numerical studies of these dif- 
ferent approaches, and a Master's degree thesis on the subject is ’n pre- 
paration. 

If we are able to conclude from this evaluation that Cliff's method 
provides a reliable test for the number of predominant modes, then we still 
must determine if there is, in fact, a relationship between this number and 
the number of shakers required to isolate modes. It is intuitively appealing 
that there should be such a relationship, but its existence has not been 
established. 

VI. Other Work in Progress 

Our studies thus far of the RFAD method have considered only hystereti cully 
damped structures. However, the results of Ibrahim and Mikulcik (1976) and 
the continuing work of Ibrahim at NASA LaRC suggest that the damping of real 
structures might be represented more accurately by a viscous model or 
possibly a combination hysteretic-viscous model. Therefore, we are extending 
the capability of our RFAD software to the case of arbitrary viscous damping. 
The theory employed to solve the forced response problen. 's that presented in 
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detail by Hurty and Rubinstein (1964) and used also by Ibrahim and Mikulcik. 

We have also begun the design of a simple experimental model patterned 
after our basic analytical model.. Fig. 1. We plan to use aluminum rods as 
the grid beams and steel rods with axes oriented vertically as the joint 
masses. To analyze this specimen accurately, we are improving our analytical 
model to account for the distributed inertia of the aluminum beams and for 
the rotational Inertias of the steel rods. Also, to design closely spaced 
modes into the specimen, we are revising the mode-merging analysis discussed 
in Section III and Appendix A to account for the constraints that exist be- 
tW'^en the translational and rotational inertias of the steel joint masses. 

We plan to run modal survey tests on this specimen using excitation 
equipment on loan to VPI&SU from NASA LaRC and subminiature accelerometers 
recently acquired by VPI&SU. We also hope to employ a Zonic/Tektronix FFT 
data processing system, two of which are being purchased by VPI&SU. 


10 



REFERENCES 


Asher, G. W. (1958). A method of normal mode excitation utilizing admittance 
measurements. Proceedings, National Specialists Meeting, IAS; Dynamics 
and Aeroelasticity (Institute of the Aeronautical Sciences), pp. 69-76. 

Asher, G. W. (1967). A note on the effective degrees of freedom of a vibrating 
structure. AIAA Journal , Vol. 5, No. 4, pp. 822-824. 

Bishop, R. E. D. and Gladwell, G. M. L. (1963). An investigation into the 
theory of resonance testing. Phil. Trans, of the Royal Society of London . 
Series A, Vol. 255, Mathematical and Physical Sciences, pp. 241-280. 

Bishop, R. E. D. and Johnson, D. C. (1960). The mechanics of vibratio n. 
Cambridge University Press. 

Cliff, E. M. (1975). Essential uncontrollability of discrete linear, time- 
invariant, dynamical systems. Proceedings of the IFAC 75, Part 1, Section 
9.3, pp. 1-6. Published by the International Federation of Automatic 
Control . 

Craig, R. R. and Su, Y. W. T. (1974). On multiple shaker resonance testing. 
AIAA Journal , Vol. 12, No. 7, pp. 924-931. 

Hurty, W. C. and Rubinstein, M. F. (1964). Dynamics of Structures . Prentice- 
Hall, Inc. 

Ibrahim, S. R. and Mikulcik, E. C. (1976). he experimental determination of 
vibration parameters from time responses. The Shock and Vibration 
Bulletin, Bulletin 46, Part 5, pp. 187-193. 

Mead, D. J. (1970). The existence of normal modes of linear systems with 
arbitrary damping. Proceedings of a Symposium on Structural Dynamics, 

Vol. I, Paper No. C5. Loughborough University of Technology, England. 

Stafford, J. F. (1976). A numerical simulation of multiple-shaker modal 
testing. M. 5. thesis, Virginia Polytechnic Institute and State 
University, Blacksburg, Virginia. 

Traill-Nash, R. W. (1958'. On the excitation of pure natural modes in aircraft 
resonance testing. J. Aero/Space Sc i., Vol. 25, pp. 775-778. 

Zarghamee, M. S. (1968), Optimum frequency of structures. AIAA Journal , 

Vol. 6, No. 4, pp. 749-750. 



APPENDIX A 


A SIMPLE PROCEDURE FOR PRODUCING 
CLOSELY SPACED VIBRATION MODES 
IN AN ANALYTICAL STRUCTURAL MODEL 


We consider here the standard n-degree-of -freedom structural dynamic 
eigenvalue problem, 

("4 [m] + W ) tr" 0 \ r = 2. n (1) 

We assume that the mass and stiffness matrices can be separated into full and 
diagonal component matrices. 


Cm] = [mf] + [v] 

[k] = [kf] + [a] 


( 2 ) 


The and matrix elements correspond to lumped physical members, which 
often are easily varied in laboratory experimental models as well as analytical 
models. 

It is necessary first to develop expressions for the derivatives of the 
natural frequencies oa^ with respect to the diagonal mass and stiffness elements, 
\i^ and aj^. The eigensolution for the rth mode of vibration satisfies Eq. 1. 
Substituting Eq. 2* into Eq. 1 and differentiating it with respect to gives 


9o) « j, 3(|) 

(-2a,^ ^ [m] - a,; [1^ 0) + (- 0 ,; [m] + [k] ) ^ = 0 (3) 

k k 

where [Ij^] is filled with zeros except for a unit element in the kth diagonal 

position. We next premultiply Eq. 3 by and use the definition of the rth 

generalized mass, M^ = [m] Finally, we simplify the equation by 
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applying the identity 

♦/ (-4 l"i) + (k) ) • 


which is simply the transpose of Eq. 1. From this procedure we obtain 




1 ^ 2 
2M^ “r ‘**kr 


( 4 ) 


By differentiating Eq. 1 with respect to Oj^ and procedi ng as above, we may 
obtain also 


^__J . 2 

2M^“r 


(5) 


The procedure leading to Eqs. 4 and 5 is due to Zarghamee (1968). 

Now, our general objective is to develop an analytical model which is 
physically reasonable with at least one pair of closely spaced modes. To do 
this, we may start with a specific numerical structural model and adjust the 
diagonal terms and Oj^, k = 1, 2, ... n. To simplify the analysis, let 
us consider adjustments only to the mass elements, \i^. The requirement of 
physical reasonableness constrains all \i^ to be positive and less than some 
maximum possible value. 

We define vector y as the n-dimensional column of \\^ elements. Then the 
difference n between any two adjacent natural frequencies may be considered a 
function of p , 

^(m) ” “j " “-j (6) 

where we choose w. ^ w. for this discussion. The Taylor series representation 
of this functional relation is 

fi(y) = ^(Wq) + (2^1, J ' (h “ + ••• 
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trtwre ts art initial state arid it is assun^d that ii- is small so 

that terins of higher order may be neglected. Hence, with = n(u) - 3 (Hq) 

and 6u = jj- 


6H = [Vft] * «u + ... 

~ Ho 

= 6|j + ... 


(7) 


where 


„ „ 3(0- 3(i)- 

g = = [ — J- 

^0 V 

itQ 


Clearly, g may be evaluated with the use of Eq. 4. 


To achieve the objective of closely spaced modes, we seek to adjust the 
masses so that 6n is negative. Mass changes 6 m must be kept small both to 
justify the neglect of higher-order terms in Eq. 7 and to maintain physically 
reasonable mass values. A simple way to maximize 6n while holding 5 m small 
is to specify 

6m = eg (8) 

where e is a dimensional number to be determined. For a given total length 
of vector 6u, the form specified in Eq. 8 clearly maximizes the inner product 
in Eq. 7. Substituting Eq. 7 into Eq. 8 gives 

i 

The change desired in n may be written as 

6n = (10) 
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where f is a negative decimal number usually ranging from 0 to -0.5. Now 


e may be determined from Eqs. 9 and 10 as 


= f 


g^g 


( 11 ) 


Hence, the procedure used to bring together two nrades of a structural 

model by adjusting its masses is as follows: 

* 

(i) Calculate the modes of the original configuration, and from 
those modes calculate with Eq. 4 the frequency-gradient matrix. 

(ii) On the basis of the original modes and the frequency-gradient 
matrix, select two adjacent modes to be brought together. 

(iii) Using the definition of Eq. 10, select a ratio f, 0> f'- -0.5, by 
which the frequency spacing is to be reduced. If this ratio is 
too large, the assumption that higher-order terms in Eq. 7 are 
negligible may be invalid, and this procedure may fail to give 
satisfactory results. 

(iv) Now calculate the mass change vector which should produce the 
desired reduction in frequency spacing. 


6m = eg 

with e given by Eq. 11. 

(v) Calculate m-j = Uq subsequently calculate the modes for 

this new configuration. The candidate modes selected in (ii) should 
now be closer together by about the amount specified in (iii). 

(vi) If it is desirable to bring the candidate modes even closer together, 
then repeat the procedure as many times as necessary, omitting step 
(ii) in all iterations after the first. 
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The theory and procedure presented here can be broadened vHthout much 
difficulty to account for such situations as specific constraints on or 
between the mass values and variations of the modes with respect to 
both lumped masses )i|^ and lumped stiffnesses 0 |^. In the latter case, 
frequency difference n would be regarded as a function of both u and a 
and both Eqs. 4 and 5 would be used. 
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_J\PPtNOIX B 

A METHOD FOR DETERMINING THE NUMBER OF MODES 
PREDOMINANT IN FORCED SINUSOIDAL RESPONSE 

Consider a structural system with n discrete degrees of freedom. We 
denote the column matrix of structural response as x and the column matrix 
of force as f. For steady-state sinusoidal forcing and response at frequency 
b), we use complex notation, 

f 

f = Re [ Fe^“^ j = F cos ut 
X = Re [ Xe^“^] 

The complex response vector X represents both amplitude and phase of response. 

In general, one can employ modal analysis with real or complex modes to 
express forced response in the series form 

X = E K (,.) i F 
r=l -r -r - 

where is the mode shape of the rth mode and K^(w) is a constant dependent 
on the generalized nass and damping of the rth mode. Hence, the complete 
n x n transfer function matrix is 

[A] = 'f'p 

r=l 

And any column of this matrix may be expressed as a linear sum of the n mode 
shape vectors, 

tr 

where (j>. is the jth element of 
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In order to simplify this discussion, let us restrict attention to 
situations for which all K and ♦ are real. Consider, for example, a system 

I ^ f 

with hysteretic damping which does not couple the undamped normal modes. 

For such a system, the real part of the admittance matrix, which represents 
in-phase or coincident response, may be expressed as 


Re[A] = [C] = Z K^(a,) 

r=i 


where 

Mr<"r - 

in which M^, and are the natural frequency, generalized mass, and modal 
structural damping, respectively, for the rth mode. With these definitions, 
each column of the coincident-response matrix can be expressed as a series sum 
of real n-dimensional column vectors, 


c, = E U“) ♦ 


-3 


r=1 


jr :r 


In the process of sinusoidal vibration testing, it is usually possible to 
measure response at p stations on the specimen, where almost invariably p < n. 
Hence, with a single shaker, say the j th , forcing at frequency w, one measures 
an incomplete column vector of the admittance matrix, 

n 

C* = K (u) (j>. (f* 

-J ' ^jr :r 

where and are p-dimens ional column vectors which lack the n-p elements of 
the corresponding C. and f . Since the 4>* are p-diniens ional , at most only p of 
the n incomplete mode shape vectors are linearly independent. It seems reasonable 
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to expect that in general any sublet of p (|)*'s will, in fact, be independent. 

At any particular frequency of excitation, we would expect some subset, 
say q, of the total modes to be predominant in the response. If this is the 
case, then each should be a linear sum of the q «ji*'s corresponding to the 
predominant modes. This suggests the following procedure for determining the 
number q: 

A. If k shakers and p motion sensors are available, measure at 
frequency w the incomplete p x k coincident-response matrix 
v^ose col umns are the 

J 


[C*] = [C* , c* . . , c*] 

B. Now analyze [C*] as follows: 

(1) Determine the single p-dimensional unit vector h-j which, 
in a sense to be defined, best represents each column 
vector as a linear sum, i.e., 

cj " j = 1, 2, . . k 

Calculate, in a sense to be defined, the error resulting 
from this representation. 

(2) Determine the pair of p-d1mensional orthogonal unit 
vectors h^ and h2 which best represent each as a 
1 inear sum, i.e., 

2 ? 

q " E C^.h., j = 1, 2, . . k 

^_1 ji-i 
2 

Calculate the error E resulting from this representation. 
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(m) Determ-fne the m p-cMmensional orthogonal unit vectors h-j , h^* 
. . h which best represent each C’J as a linear sum. 
Calculate the error e"’ resulting from this representation. 


If after q steps of this procedure we find and E^ ~ 0, 

then we may reasonably conclude the each is, with very small 
error, a sum of q orthogonal unit vectors and, therefore, that q 
modes predominate in the response at the frequency in question. 

If the procedure is carried through p steps, it is clear that 
E =0 since the p independent h. span R^. 

The theory required to carry out the calculations in B. above is next 
developed. To simplify the notation we denote the given incomplete coin- 
cident-response matrix as 

[C*] = [Z] = [z^, z^, . . ., Z(^] 

Let a basis for consist of p orthogonal unit vectors h^, h 2 , . • .. hp, 
which are unknown at this point. Hence any real p-dimensional vector can be 
expressed as the sum 


z . 





j = 1, 2, 


• 1 


k 


Since the unit h^. are orthogonal, 


7 . . 


• z . = 





2 


Hence, we may define the scalar error resulting from representing z. by only the 

3 

first m of the h^. as 
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It follows that an appropriate total error, defined in the least-squares sense, 
for all z., j = 1, 2, . . ., k, is 

M J 


-in 


r k 
= i; 
L j*i 


m \ 2 


(ep 


1/2 


‘ k km 5 

= S • Zj - E E (h. • Z.)^ 

L j=l j=l i=l ' J 


1/2 


J = i J= 

Clearly, the maximum value of the double-summation corresponds to the minimum 
E*", which is the value that we seek. It can be shown that 
k 


i?, (li • = 

J * 


where 


tT] = [Z][Z]‘ 

[T] is a p X p real, symmetric, positive-semidefinite matrix, so it has p non- 
negative real eigenvalues and p corresponding real orthogonal eigenvectors 
. It is convenient to take the to be unit vectors. It can be shown that 
the maximum values of the quadratic form involving [T] are 

ihpT]h,l -♦•[T]*. 

I- ' 'max 

Hence, in order to minimize we take h^. = , i = 1 , 2, . . . , p, and we obtain 




1/2 
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We can now reiterate and summarize the method for determining the number 
q of predominant modes at a given frequency; 

(i) Using experimental or analytical data, assemble the p x k 
incomplete coincident-response matrix [C*] based upon k 
shaker locations and p sensor locations. 

(ii) Form the p x p matrix 

[T] = [C*][C*]‘ 

(iii) Calculate the eigenvalues of [T], X., i = 1, 2, . . p. 

(iv) For m = 1, 2, . . p, calculate the error 




m 


i: 

i=l 



If there are in fact q < p predominant modes, then this should be indicated 
by the observations that and E*^ = 0. 
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